Bethe ansatz solution of zero-range process with nonuniform stationary state 
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The eigenfunctions and eigenvalues of the master-equation for zero range process with totaUy 
asymmetric dynamics on a ring are found exactly using the Bethe ansatz weighted with the station- 
ary weights of particle configurations. The Bethe ansatz applicability requires the rates of hopping 
of particles out of a site to be the q- numbers [n],. This is a generalization of the rates of hopping 
of noninteracting particles equal to the occupation immber n of a site of departure. The nonin- 
teracting case can be restored in the limit q 1. The limiting cases of the model for q = 0,oo 
correspond to the totally asymmetric exclusion process, and the drop-push model respectively. We 
analyze the partition function of the model and apply the Bethe ansatz to evaluate the generating 
function of the total distance travelled by particles at large time in the scaling limit. In case of 
non-zero interaction, g 7^ 1, the generating function has the universal scaling form specific for the 
Kardar-Parizi-Zhang universality class. 
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I. INTRODUCTION 

The Bethe ansatz U is one of the most powerful tools 
to get exact results for the systems with many interact- 
ing degrees of freedom in low dimensions. The exact 
solutions of one-dimensional quantum spin chains and 
two-dimensional vertex models are classical examples of 
its application 0]. In the last decade, the Bethe ansatz 
was shown to be useful to study one-dimensional stochas- 
tic processes 0, 01 • The first and most explored exam- 
ple is the asymmetric simple exclusion process (ASEP), 
which serves as a testing ground for many concepts of 
the nonequilibrium statistical physics j^. Yet several 
other Bethe ansatz solvable models of nonequilibrium 
processes have beeiiproposed such as the asymmetric 
diffusion models |^ [3 , generalizations of the drop-push 
model IM Hoj , and the asymmetric avalanche process 

(ASAP)ll3- 

Most of models studied by the coordinate Bethe ansatz 
have a common property. That is, a system evolves to 
the stationary state, where all the particle configurations 
occur with the same probability. This property can be 
easily understood from the structure of the Bethe eigen- 
function. Indeed, the stationary state is given by the 
groundstate of evolution operator, which is the eigenfunc- 
tion with zero eigenvalue and momentum. Such Bethe 
function does not depend on particle configuration at all 
and results in the equiprobable ensemble. Except for a 
few successful attempts to apply the Bethe ansatz to sys- 
tems with nonuniform stationary state, such as the ASEP 
with blockage or defect particle [13 , there is no much 
progress in this direction. 

On the other hand many interesting physical phe- 
nomena such as condensation in nonequilibrium sys- 
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tems Q, boundary induced phase transitions 0, 
or intermittent-continuous flow transition become 
apparent from non-trivial form of the stationary state, 
which change dramatically from one point of phase space 
to another. Typical example is the zero-range pro- 
cess (ZRP), served as a prototype of a one-dimensional 
nonequilibrium system exhibiting the condensation tran- 
sition . While its stationary measure has been studied 
in detail the full dynamical description is still ab- 
sent. 

The aim of this paper is to obtain the Bethe ansatz 
solution of the ZRP. The paper is organized as follows. 
In section II we use the Bethe ansatz to solve the eigen- 
function and eigenvalue problem for the master equation 
of the ZRP and show that its applicability requires the 
rates of hopping of particles out of a site to be the q- 
numbers. We show that with this choice of the rates the 
model is equivalent to the q-boson totally asymmetric 
diffusion model. In section III we obtain the partition 
function of the ZRP with the rates obtained and evalu- 
ate some stationary correlations. In section IV we apply 
the equations obtained from the Bethe ansatz solution 
to derive the expression for the generating function of 
the total distance travelled by particles in the large time 
limit. We summarize the results in section V. 



II. MASTER EQUATION OF ZERO RANGE 
PROCESS 



Let us consider the system of p particles on a ring con- 
sisting of N sites. Every site can hold an integer number 
of particles. Every moment of time, one particle can leave 
any occupied site, hopping to the next site clockwise. 
The rate of hopping u{ni) depends only on the occupa- 
tion number Ui of the site of departure i. The stationary 
measure of such a process has been found to be a product 
measure 0|, i.e. the probability of configuration {n;}. 
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specified by the occupation numbers {rii, n2, . . . , Hn}, is, if n > 0, and /(O) = 1. 
up to tlie normalization factor, given by the weight 



N 



W{{n,}) = Y[fin,), 



(1) 



where one-site weight /(n) is 



n 

/H ^ n ^ 

^ ^ mm) 

m— 1 ^ ^ 



(2) 



Let us consider the probabihty Pt{ni, . . . ^un) for N 
sites to have occupation numbers ni, . . . , njv at time t. 
It obeys the master equation defined by the dynamics 
described. 



dtPt{ni, . . . .un) = ^ + 1)P4(. . . ,nfe_i + l,nfe - 1, . . .) - M(nfc)Pt(rii, . . . ,nAr)) (3) 

k = 1 



Here, cyclic boundary condition, n_i = mat, is imposed. 

Another way to specify the configuration of sys- 
tem is to use the set of coordinates of p particles 
{xi\ — {xi, . . . ,Xp\, rather than occupation numbers 
{ni, . . . , riAr}, the two ways being completely equivalent. 
While the explicit form of the master equation in such 
notations is more complicated, it turns out more ap- 
propriate for an analytic solution. The main idea of 
the solution is to use the Bethe ansatz for the function 
P^{xi,...,Xp) related with the solution of the master 
equation Pt{xi, . . . , Xp) as follows 

Pt{xi,...,Xp)^W{{n,})P°{ )• (4) 

A. Two-particle sector 

To explain the details we first consider the case p = 2 
subsequently generalizing it to the case of arbitrary p. 
Without lost of generality we can define 

u(l) = l, u(2) = u>0. (5) 

Now, we are going to show that the solution, 
Pf{xi, X2), of the master equation for noninteracting par- 
ticles, 

dtPt^{xi,X2)^ Pt°{xi-l,X2) + P°(X1,X2-1) 

- 2P,°ixi,X2), (6) 

is related with the solution Pt{xi, X2) of the master equa- 
tion for the ZRP in the domain xi < X2 through the re- 
lation (0J, provided that the former satisfy certain con- 
straint. Indeed, when {x2 — Xi) > 2, the equation for 
probability Pt{xi,X2) for the ZRP coincides with one for 
noninteracting particles ®. For xi — X2 — 1 — x the 
equation for the ZRP 

dtPt{x,x+l) ^ Pt{x-l,x)+uPt{x,x)-2Pt{x,x+l) (7) 



can be obtained from (jSj) if we define 

Pt{x,x)^-P°ix,x), (8) 
u 

which is consistent with I^J. In the case X2 — xi — x 
the equation ((HJ contains the term P^{x,x — 1) which 
is beyond the "allowed" region xi < X2 and, thus, does 
not carry any physical content. To restore the correct 
equation for the ZRP 

dtPt{x, x) = Pt{x -l,x)- uPt{x, x), (9) 

we can redefine this term to compensate the difference 
between the equations (O and As a result we get 
the following constraint on P^{xi, X2) 

P^{x,x-1) = iu~l)Pt"ix-l,x)-{u - 2) P,"ix,x). (10) 

Thus, the solution of the free equation ©, which satis- 
fies the constraint (|10|) . gives the solution of the master 
equation for the ZRP in the domain xi < X2- Now, we 
can use the Bethe ansatz for the eigenfunction of the free 
equation (0 

xa) = e^* (Ai,2zr"^^2""' + ^2,i^r"^^2""0 ■ (11) 

Its substitution to the equation ^ results in the expres- 
sion for the eigenvalue. 

A = zi + Z2 - 2. (12) 

The ansatz (|11|) to be consistent with the constraint H1U|I , 
the amplitudes Ai^2 and A2.1 should satisfy the relation 

^1,2 _ (2 - m) - (1 - U)Z2 - Zl ^^^^ 

^2,1 (2 - U) - (1 - U)Z1 - Z2 ' 

which together with the cyclic boundary conditions, 
P^{xi,X2) = Pf{x2,xi + N), results in the system of 
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two algebraic equations. The first one is the following 

-N ^ {2-u)- {l-u)z2- Zi 

while the second can be obtained by the change zi <-> Z2. 

B. Many- particle sector 

To generahze these resuhs to the case of arbitrary p 
let us consider the configuration with two neighboring 



sites (s — 1) and x having occupation numbers m and n 
respectively. Let us explicitly write down the terms of 
the master equation corresponding to the transition into 
and from this configuration due to a particle jump into 
and from the site x respectively, 



dtPt{. .-Ax- ir, (x)", ...) = ...+ + i)Pt(. . . , (x ~ iT^\ {xr-\. . .) 

-«(n)Pt(...,(.T-l)™,(a:)",... ) + .... (15) 

Here (a;)" denotes n successive arguments x of the function ^4(2^17 ■ ■ ■ , Xy\ In terms of P^(x\^ . . . , Xp), which is related 
to Pt{x\^ . . . , Xp) according to this equation looks as follows 

9tP0(. . . , (x - 1)", (x)", ...) = ... + uin) X [P0(. . . , (x - l)™+\ (x)"-\ . . .) 

-P°(...,(x-l)'",(x)",... )]+... (16) 

Note, that in this form the coefficient before the term in square brackets is equal to u(n), i.e. does not depend on the 
number to of particles in the site (x — 1). Thus, the r.h.s of the master equation expressed through P(*^(xi, . . . ,Xp) 
is the sum of one-site factors similar to those in (|16|l for all nonempty sites. Equating such term with one from the 
equation for noninteracting particles, 

V 

dtP^ixi, . . . , xp) = ^ (P°(. . . ,x, - 1, . . .) - P0(. . . ,x„ . . .)) (17) 
we obtain the following constraint for P°(xi, . . . , Xp), 

n 

{u{n) 1)P0(. . . , (x), (x + ir-\ ..:)-Y^ P°(. . . , (x + iy-\ (x), (x + 1)"-^ . . .) 

-(K(n)-n)pO(...,(x + l)",...) = 0. (18) 
I 



Such relations for all n are to be understood as a redefini- 
tion for the terms outside of allowed region xi < . . . < Xp. 
The Bethe ansatz to be applicable, the relation H18|) 
should be reducible to the two particle constraint (|10|l . 
This can be proved by induction. To this end, we assume 
that similar relations including u(k) are reducible for all 
k < n. Then starting from the relation (|18|l . which in- 
cludes rate u{n), we apply (|l()|l to every pair (x,x — 1) 
of arguments of the function P^{. . .) under the sum and 
require the result to be similar relation for u{n— 1). We 
obtain the following recurrent formula for the rates 

w(n) = 1 - (1 - u)u(n- 1), (19) 

which can be solved in terms of g-numbers 

u{n) = = , (20) 



where 

q = u-l. (21) 

Further generalization of two-particle results is straight- 
forward. We use the Bethe ansatz for the eigenfunction 
PjP(xi, . . . , Xp) of the equation ((T7|l 

P°(xi,...,xp) = e^* ^ A{,,_,^}\{z-:^ (22) 

{cri,...,(Tp} 1=1 

Here zi, . . . ,Zp are complex numbers, the summation is 
taken over all p\ permutations {cti, . . . , (t„} of the num- 
bers (1, . . . ,p), and the coordinates of particles are or- 
dered in the increasing order xi < X2 < . . . < Xp. 
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Substituting ^ into (UTj) we obtain the expression 
for the eigenvalue, 



(23) 



The numbers (zi, . . . , Zp) are to be defined from the Bethe 
equations. 



-N 



(2 — 7i) — (1 — u)Zj — Zi 

(2 - m) - (1 - u)zi - Zj 



(24) 



which follow from the condition of compatibility of cyclic 
boundary conditions with the constraint H10|l . 

For the sake of convenience in the following discussion 
we use the parameter q defined in H21|l rather than u. 
We should note that the appearance of q-numbers as the 
condition of the Bethe ansatz integrability is not unex- 
pected. The notion of q-deformation naturally appears in 
context of Bethe ansatz solvable models characterized by 
the trigonometric R-matrix. In algebraic language this is 
the consequence of the of the fundamental Yang-Baxter 
equation which leads to q-commutation relations between 
the local operators constituting the transfer matrices |0| . 
One of such models, q-boson totally asymmetric diffusion 
model 0, turns out to be directly related to the model 
we consider. In order to see the correspondence, let us 
formally write the distribution Pt (C) as a vector of state 



|P*(C))=^P,(C)|C), 

{C} 



(25) 



where C is a configuration of particles on the lattice, 
C = {ni, . . . , n^r}, and the summation is over all configu- 
rations. Consider the algebra generated by the operators 
Bj , B^ , N'j , which act on the occupation number nj > 
of each site j of the lattice as follows: 



Bj \nj) = \nj - 1) 



B+\n,) = [nj + l]q\nj + l) 



(26) 
(27) 
(28) 



Afj \nj) = Uj \nj) , 

The state |0) plays the role of vacuum state 

Bj \0j) = 0. (29) 

Then, the master equation Q with the rates given by 
(|20|l can be written in the from of the imaginary time 
Schrodinger equation 



dt\PtiC))^~Ji\PtiC)) 



(30) 



where the hamiltonian H is given in terms of the opera- 
tors (??iinii 



ll = -J2iB;_,B,-B+B, 



(31) 



One can directly check that the operators Bj , Bj' , A/'j sat- 
isfy the following commutation relations 

[Bj,B+] - q^^S,,, (32) 
[JVk,Bj] = -BjSjk (33) 
[Mk,B+] = B+6jk- (34) 

These commutation relations and the hamiltonian l|31|) 
give us, up to the change of notations q — > (?~^, the defi- 
nition of the q-boson totally asymmetric diffusion model. 
Obviously, the dynamical rules of the ZRP with the rates 
given by H2U|) is nothing but the explicit realization of 
this hamiltonian. Its integrability has been shown in 
via algebraic Bethe ansatz and two particle diffusion on 
the infinite lattice has been considered. Note that while 
in the q-boson totally asymmetric diffusion model 
was initially defined in terms of q-boson operators, we 
started from the ZRP with arbitrary rates and come to 
q-numbers as the integrability condition. 

Let us first take a qualitative look at the behaviour of 
the ZRP resulted by the choice (|20() of the rates u{n) for 
different values of q. In the limit q 1 the g-numbers 
degenerate into simple numbers. Therefore, the rates are 
given by u{n) = n, which corresponds to the diffusion 
of noninteracting particles. The Bethe equations in this 
case decouple to the form zj^ = 1 as is expected in non- 
interacting case. In the domain q > 1 the rates u{n) 
grow exponentially with n, which corresponds to the in- 
teraction between particles effectively accelerating free 
diffusive motion, i.e. the higher density of particles the 
faster is their mean velocity. In the limit q ~* oo the 
model is equivalent to n = I drop push model ^ , which 
is confirmed by the same form of Bethe equations ^ . In 
the domain < q < 1 the rates u{n) grow monotonously 
from {q-\-l) for n — 2 to 1/(1 — g) for n oo, resulting in 
the interaction between particles slowing down the parti- 
cle flow comparing to the one of noninteracting diffusing 
particles. When q = 0, all the rates do not depend on 
the number of particles, i.e. u(n) = 1. This case, (also 
referred to as phase model 113), can be mapped on the 
totally asymmetric ASEP by insertion of one extra bond 
before every particle. At last, in the domain — I < g < 
the rates u{n) also saturate to the constant 1/(1 — q) 
with growth of n, though oscillating around this value. 
As it has been mentioned above, the ZRP served as an 
example of the nonequilibrium system with the conden- 
sation transition. In our case however the condensation 
is absent as the rates u{n) defined above do not satisfy 
Evans criteria, according to which the condensation in 
the ZRP occurs if the rates saturate to a constant slower 
than 2/n. 

It is interesting, that the recursion relation (|f 9|) rewrit- 
ten in terms of the quantity (1 — u{n)) coincide with one 
for the toppling probabilities /i„ , imposed by the require- 
ment of the Bethe ansatz integrability in the ASAP. In 
fact, the ASAP can be represented as a special case of 
the discrete time ZRP viewed from the reference frame 
moving together with an avalanche. This situation, how- 
ever, is quite different from one considered here. In the 
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moving reference frame all particles hop definitely to the 
next site except of one from an active site, which is the 
only site with multiple occupation. In that case quan- 
tity 1 — /i„ plays the role of probability of hopping of 
a particle out of this site. This dynamics leads to the 
picture similar to the ZRP on an inhomogeneous lattice 
with one attractive site. Such ZRP was shown to exhibit 
the condensation transition. In terms of the avalanche 
processes it is the transition from the intermittent to con- 
tinuous flow. Despite the nonuniform stationary state of 
the ASAP in discrete time j its Bethe ansatz solution 
was based on the continuous time picture considered on 
the ensemble of equiprobable configurations with at most 
one particle occupation. The avalanches were accounted 
in the rates of transitions between these configurations, 
generating infinite series in the master equation. 

To use the eigenfunctions obtained for the construction 
of particular solutions one should first question if they 
form complete orthogonal basis. In general this ques- 
tion is not easy to answer, as the set of solutions of the 
Bethe equations is not known. Some arguments have 
been given, for the the totally asymmetric exclu- 

sion process due to its connection with the asymmetric 
six- vertex model[2^. The long time characteristics of the 
process can, however, be analysed without discussing this 
question. To this end, we can use the properties of the 
largest eigenvalue for slightly modified equation, which 
describes the generating function of total distance trav- 
elled by particles, l2^. The advantages of this approach 
are first, that the uniqueness of the largest eigenvalue is 
guarantied by Perron-Frobenius theorem. Second, corre- 
sponding solution of the Bethe equations can be easily 
identified as it corresponds to the stationary state of the 
process. To give an example of application of the above 
results we perform this analysis in the section IV. 



for \q\ < 1 



Fiz) 



E 

»i=0 



(^(l-g))" 



= (2(1 - q);q)c 



(37) 



and for \q\ > 1 



F(z) 



t^'^'lTC'^-"" (38) 



i<q- 



{q-';q-')n 
l);q~')oo- 



Here, we used the notation {a;q)n = nfe=o(l ~ '^9'^) fo'^ 
shifted g-factorial. The above g-series are known as q- 
analogs of the exponential function e^, which can be re- 
stored in the limit g — > 1. The presence of g-analogs 
of the functions, which appear in the case of noninter- 
acting particles, is a direct consequence of the replace- 
ment of simple numbers by g-numbers in the expression 
of the rates of hopping. In the thermodynamic limit, 
N 00, p 00, p/N — p, the integral H36() can be cal- 
culated in the saddle point approximation. The equation 
for the saddle point zq. 



P = zo^og' F{zo), 



(39) 



contains the logarithmic derivative of F{z), which can be 
evaluated using the product form of F{z) H37I38(I . As a 
result we obtain the following relations between zq, p and 

q- 



P = zoi^ - q) 



^ 1 

n=0 



2o(i - q)q'^- 



(40) 



III. STRUCTURE OF THE STATIONARY 
STATE 

Before going to the results of the analysis of the Bethe 
equations, let us first look at the structure of the sta- 
tionary measure of the model. It is characterized by the 
partition function 



00 

Z{N,p) = V (5(711, 



N 



■ + nN-p)Y[fin^), (35) 

1=1 

where one site weight /(n) = l/[n]q!, defined in is 



{jii 



,71 



, = 1} 



expressed through the g-factorial 



n 



I — 



the limit q —> 1, 5- factorial turns into simple factorial, as 
it should be in the noninteracting case. The sum in H35|) 
can be presented in form of the contour integral 



ZiN,p) 



1 

2TTi 



{F{z)) 



N 



-dz. 



(36) 



where the series F{z) = X^J^Lo /('^)2" can be summed to 
the infinite product due to the g-binomial theorem [24| : 



^ I - zo{q 1 - I)g " 

for |(?| < 1 and |(7| > 1 respectively. Below the same equa- 
tions will appear in a different context from the analysis 
of Bethe ansatz equations. Then, the partition function. 



N 



(42) 



can be used to calculate stationary state correlations such 
as the speed, v = (ii(n)), i.e. the average hopping rate 
out of a site 



Z{N,p-l) 



(43) 



or the probability distribution of the number of particles 
in a site 



[nW-FizoY 
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IV. THE LONG TIME BEHAVIOUR FROM 
THE BETHE EQUATIONS 

To obtain any results beyond the stationary correla- 
tions one needs to analyze the above Bethe ansatz solu- 
tion. Since similar analysis has been done several times 
before [U IIE HE 113 j we do not give the detailed calcu- 
lations here. Instead we outline the main points of the 
solution to emphasize the connection with the formulas 
obtained from the analysis of stationary measure. 

Consider the generating function ^t{C) = 
T.Y=oPt{C,Y)e^^ , where Pt{C,Y) is the joint prob- 
ability for the system to be in the configuration C at 
time t the total distance travelled by particles being 
Y . The only difference of the equation for Tt{C) from 
the equation |(2Jl is the coefficient before the first 
term under the sum in the r.h.s., which corresponds to 
the increasing of travelled distance by one due to the 
hopping of one particle. At large time, t ^ oo, the 
behaviour of the generating function of the distance Yt 
travelled by particles up to time t, {e"^^*) = J2c ^ti^)' 
is determined by the largest eigenvalue of the equation 
for J^tiC). 

Ao(7)= lim i^i^ (45) 

Using the ansatz (|4I22(I for the eigenfunction of the equa- 
tion for Tt{C) we can repeat all the above arguments. 
Then, if we make the variable change Xi = 1 — Zi, the 
eigenvalue and the Bethe equations will simplify to the 
following form. 

A(7) = -Ea;. (46) 

i=l 

e^^(l - rrO""^ = (-l)P-' n ^^^^ (47) 

In these variables the r.h.s. of the Bethe equations co- 
incides with one for the partially asymmetric ASEP and 
the ASAP. This allows us to modify the techniques de- 
veloped for the analysis of these processes. 

In the thermodynamic limit, iV — s- oo,p — > oo, p/N = 
p, we assume that the roots of the Bethe equations (|47f) 
are distributed in the complex plain along some contin- 
uous contour F with the analytical density R{x), so that 
the sum of values of a function f(a:) at the roots is given 

by 

P r 
Vf(x,) i{x)R{x)dx. (48) 

^=l -^r 

After taking the logarithm and replacing the sum by the 
integral, the system of equation H47|) can be reformu- 
lated in terms of single integral equation for the den- 
sity. The particular solution corresponding to the largest 



eigenvalue is specified by the appropriate choice of branch 
of the logarithm. Then the integral equation should be 
solved for a particular form of the contour, which is not 
known a priori, and being first assumed should be self- 
consistently checked after the solution has been obtained. 
In practice, however, analytical solution is possible in the 
very limited number of cases. Particularly, one, corre- 
sponding to the contour closed around zero, yields the 
density 

for \q\ < 1 and 

for l^l > 1. This case corresponds to 7 = and 
hence A(7) = 0. Since R^^\x) is analytic in the ring 
< < 1, the integration of it along any contour closed 
in this ring does not depend on its form. Therefore, to fix 
the form of the contour additional constraints are neces- 
sary. Such constraint appears if we require that the den- 
sity preserves its analyticity when 7 deviates from zero 
and the contour becomes discontinuous. This constraint 
implies that the density R^^'>{x) vanishes at the break 
point Xc, which is a crossing point of the contour F and 
the real axis. 

R^°Hxc) = (51) 

This equation was first obtained by Bukman and Shore as 
a conic al p oint condition for the asymmetric six-vertex 
model [23. It is remarkable, that after the resumma- 
tion of series in H49I5U|) , the equation (|51|l coincides with 
eas. (|41l40f) up to the replacements 

Xc = 2:0(1 - (?) and Xc ^ zo{q^'^ - I) (52) 

for l^l < 1 and \q\ > 1 respectively. Remember that zq 
was shown to coincide with the speed v. It will be clear 
below that the same relation between v — {X'q\^=o)/N 
and Xc follows directly from the expression for Ao(7) ob- 
tained from the Bethe equations, without appealing to 
the partition function. 

The further analysis is related with the calculation of 
finite size corrections to the above expression of {x) , 
which makes possible to probe into the nonzero values 
of 7. This can be done with the help of method devel- 
oped in 1^ [2^. Its essential part, is the construction 
of the inverse expansion Z~^{x) to the function of the 
number of roots Z{x) = R{x)dx near the break point 
Xc- Since the derivative of Z{x) in the thermodynamic 
limit vanishes (see ea. (|51|) '). its inverse expansion reveals 
the square root singularity, which in its turn becomes 
the origin of terms in the finite size expansion of 

R{x). As a result we obtain the following parametric 
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dependence of R {x) on 7, both being represented as 
functions of the same parameter C. 



Rs = 



1 1 



N3/2 2Tri 1 - q\s\ 



^(2N —^i — V2r^'"+^^^^ 



n=0 
1 



(53) 



Ar3/2 

00 



/ « V r(n + 1) c2„+i , . , - 

) — —1 ^Lwi(C), (54) 



Here i?^ and /^i"'' are the Loraunt coefficients of R{x) 
and i?'-''-' (x) respectively defined as follows 



R{x)= J2 



,s+l 



(55) 



C2n+i s and C2n+i are the coefficients of x" in 
Lo o-kX^) and log (X^felo "fe^:'') respectively, where 
a„ are the coefficients of the inverse expansion Z~^{x) 
near the point Xc- The location of Xc is to be self- 
consistently defined from the equation R {xc) — 0. For 
the first three orders of 1/ ^/N expansion the coefficients 
a„ can be obtained from the inverse expansion of zero 
order function Z'^°\x) = f R^^\x)dx, while R^^^x) has 
been obtained above. To evaluate the sum over the roots, 
one needs to integrate along the contour F, which can be 
obtained from the initial closed contour by cutting out 
small segment connecting two roots closest to the point 
Xc- For a function which is defined by the expansion 



f(x)=^f. 



(56) 



s=l 



this yields 



^f(x,) = 27rziV^q±'^f,i?, 



(57) 



minus and plus in power of g being for \q\ < 1 and \q\ > 1 
respectively. Finally the point Xc enters to all the results 
trough the coefficients C2n+i,s and C2n+i- It is related 
with the physical quantities through eas. (|40l41l43l52|l . 
We use this relation, to write the final results as a func- 
tion of speed v, density p, and the parameter q. Below 
we give the expression for the largest eigenvalue in the 
scaling limit jN^^^ = const, N ^ 00 

Ao(7) -iVw7 + fciG(fc27). (58) 
Here the function G(x) has the following parametric form 



G{x) 

X 



-Li5/2(C) 

-Li3/2(C^) 



(59) 
(60) 



with LiQ(a;) = J^kLi 1^°^ is the function of polyloga- 
rithm and the constants fci , k2 are 



1 



vl{l- q) <(t'(l-9)) 



iV3/2 V Stt [g'q{v{l - q))fl^ 
k2 = N^/^^27Tvil~q)g'^ivil-q)) 

for \q\ < 1, and 

1 Iv/il-q-^) 9i/JHl'' -^)) 



Ar3/2 



Stt 



5/2 



(61) 
(62) 

(63) 



k2 = -N^'^^2^v{l - q-^)g'y^{v{q-^ - 1)) (64) 
for \q\ > 1, were 

5^(-) = Et^- (65) 

ri=l 

The scaling form of the function G{x) was suggested to 
be universal for Kardar-Parisi-Zhang (KPZ) universality 
class '28, "29]. Using the generating function obtained one 
can evaluate all the cumulants of the travelled distance 



lim 

t — >oo 



5"Ao(7) 



(66) 



7=0 



The large deviation function, ldf(a;) — 
\imt-too ^og P{Yt ^ x) /t, can be also obtained as a 
Legendre transformation of Ao(7). 



V. SUMMARY AND DISCUSSION 

To summarize, we apply the Bethe ansatz to solve zero 
range process with the totally asymmetric dynamics on 
a ring. The cigcnfunctions of the master equation have 
the form of the Bethe function weighted with the sta- 
tionary weights of corresponding particle configurations. 
The requirement of Bethe ansatz integrability leads to 
the special choice of the rates of hopping of particle out 
of a site. The rates should be g-numbers, [n]g, generaliz- 
ing the case of noninteracting diffusing particles, where 
the rate is equal to n, the number of particles at a site of 
departure. The noninteracting case can be restored in the 
limit g 1. Two other limiting cases, q = and g — > 00 
reproduce well known totally asymmetric exclusion pro- 
cess and drop-push model respectively. The case of gen- 
eral q is shown to be equivalent to the q-boson totally 
asymmetric diffusion model. Continuing analogy with 
noninteracting case, we show that many quantities char- 
acterizing the stationary state correlations of the model 
turns out g-analogs of corresponding functions appearing 
in the noninteracting case. To provide an example of ap- 
plication of the Bethe ansatz solution obtained, we derive 
the expression for the large time limit of the generating 
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function of cumulants of the total distance travelled by 
particles. It has a universal form specific for KPZ uni- 
versality class. The question whether the q-boson totally 
asymmetric diffusion model belongs to KPZ class was ad- 
dressed in concluding remarks in j7|]. The result (|58I60() 
is an argument in favour of this assumption. 

In connection with above results the following ques- 
tions appear. First, is it possible to generalize the pro- 
posed combination of the Bethe ansatz with stationary 
weights to any other processes with nonuniform station- 
ary state, say asymmetric exclusion process with paral- 
lel update? The consideration of the associated vertex 
models is also of interest. The different weights of vertex 
configurations depending on the order of vertices would 
result in the appearance of non-local interaction. Sec- 
ond, can one apply the matrix product method to study 
the exclusion process with long range interaction asso- 
ciated with zero range process considered here to probe 
into spatially inhomogeneous situation, e.g. at the open 
chain. Appearance of q-numbers seems to be an indica- 



tion of this possibility. We expect that the matrix prod- 
uct ansatz should be again appropriately weighted with 
stationary weights of some homogeneous system. The 
consideration of such a system is attractive, as the extra 
parameter q could result in a reacher phase diagram com- 
paring to the usual totally asymmetric exclusion process. 
Third, it is interesting to establish correspondence of the 
large scale behaviour of the proposed process with KPZ 
equation. 



Acknowledgments 

The author is grateful to V.B. Priezzhev and V.P. 
Spiridonov for stimulating discussion. The work was sup- 
ported partially by FCT grant SFR/BPD/11636/2002 
and by the grant of Russian Foundation for Basic Re- 
search No.03-01-00780. 



H. Bethe, Zeitschrift fur Physik 71 205 b (1931) [16 
R.J. Baxter 1982 Exactly Solved Models in Statistical [17 
Mechanics Academic Press 

L.H. Gwa and H. Spohn Phys. Rev. A46 844 (1992) [18 

G.M. Schiitz Phase Transitions and Critical Phenomena 

19 ed Domb C and Lebowitz J (London: Academic) [19 

(2001) 

B. Derrida Phys. Rep. 301, 65 (1998) [20 
T. Sasamoto and M. Wadati Phys. Rev. E 58 (1998) 

4181 [21 
T. Sasamoto and M. Wadati J. Phys. A: Math. Gen. 31 [22 
(1998) 6057 [2i 
G.M. Schiitz, Ramaswamy R and Barma M J. Phys. A; 
Math. Gen. 29 837 (1996) [24 
M. Alimohammadi, V. Karimipour and M. Khorrami 
Phys. Rev. E 57 6370 (1998) 
M. Alimohammadi, V. Karimipour and M. Khorrami J. [25 
Stat. Phys. 97 373 (1999) [26 
V.B. Priezzhev, E.V. Ivashkevich, A.M. Povolotsky and [27 

C. -K. Hu, Phys. Rev. Lett. 87 084301 (2001) 

G. Schiitz J. Stat. Phys. 71 471 (1993) [28 
B. Derrida and M.R. Evans J. Phys. A; Math. Gen. 32 
4833 (1999 ) [29 
M.R. Evans Brazilian Journal of Physics 30 42 (2000) 
J. Krug Phys. Rev. Lett. 67 1882 (1991) 



G. Schiitz and E. Domany J. Stat. Phys. 72, 277 (1993) 

A. M. Povolotsky, V.B. Priezzhev and C.-K. Hu Phys. 
Rev. Lett. 91, 255701 (2003) 

S. Grosskinsky, G.M. Schiitz, H. Spohn J. Stat. Phys. 
113 389 (2003) 

L.D. Faddeev, Les-Houches lectures, preprint 
hep-th/9605187 (1996) 

N.M. Bogoliubov, A.G. Izergin and N.A. Kitanin, 

Nucl.Phys B 516 501 (1998) 

D. Kim Phys. Rev. E 52 3512 (1995) 

R. J. Baxter J. Stat. Phys. 108 1 (2002) 

B. Derrida, J. L. Lebowitz, Phys. Rev. Lett. 80, 209-213 
(1998) 

G.E. Andrews, R. Askey and R. Roy, Special Functions, 
Encyclopedia of Mathematics and Its Applicatons (Cam- 
bridge University Press, Cambridge, 1999) 
D.J.Bukman and J.D. Shore, J. Stat. Phys. 78 1277 (1995) 
D.-S. Lee and D. Kim, Phys. Rev. E 59, 6476 (1999) 

A. M. Povolotsky, V. B. Priezzhev and Chin-Kun Hu, J. 
Stat. Phys. Ill 1149 (2003) 

M. Kardar, G. Parisi and Y.-C. Zhang Phys. Rev. Lett. 
56 889 (1986) 

B. Derrida and C. Appert, J. Stat. Phys. 94 1 (1999) 



